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Abstract 

We investigate scenarios in which dark matter interacts with the Standard Model 
primarily through electroweak gauge bosons. We employ an effective field theory frame- 
work wherein the Standard Model and the dark matter particle are the only light states 
in order to derive model-independent bounds. Bounds on such interactions are derived 
from dark matter production by weak boson fusion at the LHC, indirect detection 
searches for the products of dark matter annihilation and from the measured invisible 
width of the Z°. We find that limits on the UV scale, A, reach weak scale values for 
most operators and values of the dark matter mass, thus probing the most natural 
scenarios in the WIMP dark matter paradigm. Our bounds suggest that light dark 
matter (m x <> mz/^ or m x 100 — 200 GeV, depending on the operator) cannot 
interact only with the electroweak gauge bosons of the Standard Model, but rather 
requires additional operator contributions or dark sector structure to avoid overdosing 
the universe. 
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1 Introduction 



Evidence for a substantial particle dark matter component in our galaxy is by now quite con- 
vincing. The strength and nature of non-gravitational DM interactions with the Standard 
Model (DM-SM interactions) are unknown, but sensitivities of many currently operating 
experiments are at the level of predicted signals for many well-studied DM models. Possi- 
bilities range in principle from just beyond current experimental bounds to levels that are 
impossible to probe directly, although the predominant weakly-interacting-massive-particle 
(WIMP) dark matter paradigm naturally provides for the observed DM relic abundance, 
while predicting DM-SM interactions mediated via couplings of approximately weak interac- 
tion strength to mediators with approximately weak scale masses. As is well-known, WIMP 
DM may be observed in three possible ways: (i) directly, through it's interactions with nuclei 
(and/or electrons) in underground detectors, (ii) indirectly, through it's self-annihilation into 
Standard Model (SM) particles in space and, lastly, (in) through it's production at colliders 
such as the LHC, appearing as an excess of missing transverse energy (MET). 

The study of WIMP DM has been historically dominated by "top-down" studies based 
on new physics frameworks (e.g. supersymmetry, extra dimensions, etc.) that are primarily 
meant to solve the gauge hierarchy problem. Such theories naturally posit new particles with 
weak-interactions and masses not far from the weak scale and more or less automatically 
provide WIMP DM candidates. In the current era much effort (rightly) continues to be 
devoted to such "top-down" studies (e.g., [lj-[6j), wherein the variety of phenomenology 
available in these frameworks can be studied in the context of fully-developed UV complete 
theories. 

In quite an orthogonal direction a more "bottom-up" approach to studying DM- 
SM interactions has recently taken shape [7j-[T7j. In the case where the DM particle is 
lighter than the degrees of freedom which mediate DM-SM interactions one can describe 
such interactions in an effective field theory (EFT) framework. In this picture the SM 
particles and the DM particle are the only light degrees of freedom in the theory and DM-SM 
interactions are described in terms of contact operators. At any fixed naive scaling dimension 
there are a limited number of such operators that respect Lorentz and gauge invariance so 
that, under the assumption that one particular operator in such a set is the dominant 
interaction channel, one can systematically derive bounds in a fairly model-independent 
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fashion. This approach lends itself particularly nicely to studies of the complementarity of 
different classes of DM search experiments since, in the contact operator approximation, the 
very same coupling that dictates DM elastic scattering is, to some extent, also that which 
dictates DM pair production at colliders and that which determines the energetic products 
of DM annihilation in astrophysical dark matter halo^J This approach has been used to 
great effect in the works [7j[8j[9j[10j[llJ, in which it was shown that if the EFT description 
holds, the irreducible collider signals provided by DM pair production with associated initial 
state radiation (mono jets) can place bounds on DM interactions with colored SM particles 
that are competitive with dark matter direct detection experiments, especially for light 
(wdm ^ lOGeV) dark matter. This approach has since been developed in many directions, 
e.g., to include more general classes of operators [12] [13] [2] and to include a wider variety 
of collider and DM searches |I5] [K] [IT] [E] . 

In this work we investigate the interactions between dark matter and the electroweak 
gauge bosons of the Standard Model. One reason to study this class of operators is simply 
that, ignorant as we are about DM-SM couplings, the most visible interactions (DM inter- 
actions with colored SM particles) may be suppressed relative to DM interactions with SM 
vector bosons. Besides this, such interactions may also indicate the extent to which DM 
is related to electroweak symmetry breaking and the extent to which the WIMP paradigm 
hold^J In "top-down" models the WIMP DM candidates are typically new partners of the 
electroweak gauge bosons (such as the gauginos in SUSY and the KK photon of univer- 
sal extra-dimensions models) or even of the higgs sector particles themselves (higgsinos in 
SUSY). As such this class of operators may provide information that is highly complementary 
to studies of WIMP dark matter in the context of UV complete theories of dark matter. 

To begin we will describe our effective operator description of DM interactions with 
electroweak gauge bosons in detail (Section [2]), focusing on fermionic DM and considering op- 
erators of naive scaling dimension d < 7. To set bounds on DM interactions with electroweak 
gauge bosons we consider weak boson fusion (WBF) searches for production of dark matter in 

3 In practice there are several caveats to this logic. For example, while the operator description may fairly 
generically apply to non-relativistic DM scattering and annihilation it may not be appropriate for typical 
LHC events with <^ 1 TeV. It may also be the case that different operators are dominant in the different 
classes of experiments (e.g., a DM-quark interaction may dominate elastic scattering while a DM-lepton 
interaction may provide the dominant DM annihilation channel). 

4 Of course the non- observation of DM-gauge boson interactions cannot negate the WIMP picture by 
itself, as the "WIMP miracle" is actually in effect in far more diverse scenarios [19] [20] . 
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8 TeV and 14 TeV LHC searches (Section |3|). We then consider DM indirect-detection bounds 
that can be obtained from current observations of dwarf spheroidal galaxies, monochromatic 
7-ray searches and from measurements of the cosmic-ray antiproton spectrum (Section [3]). 
Finally we combine collider and indirect-detection limits and discuss our results. As will be 
seen below the combined reach of these searches is quite significant, ~ 1 TeV, for a wide 
range of \ masses. The discussion presented here is related to the works [14J [21J [22j [23j [24J. 

2 Dark Matter Effective Theory 

In this section we first present a general discussion of the philosophy and structure of our 
EFT model of DM interactions with elect roweak gauge bosons, which leads us to a list of 
operators that will be the main focus of this paper (Table [I]) and their detailed description. 

We assume that the relevant light degrees of freedom include only the usual matter 
content of the SM and a dark matter particle x- Here we consider only Dirac or Majorana 
fermionicj^] x, though generalization to scalar x' s would be a straightforward extension of 
this analysis. We consider operators of naive scaling dimension d < 7 coupling x to the 
SM bosons 7, Z° and W ± (or equivalently B° and W a ), which we may generically refer to 
as V. In general we have contact operators not only of the 4-point xxYV topology, but 
also operators which generate 3-point xxY couplings to neutral gauge bosons, the latter of 
which (as we will discuss at length in the following sections) lead to substantial alterations to 
our phenomenological expectations. We require all operators to satisfy U(1)em invariance, 
but not necessarily invariance under the unified electroweak SU(2) L ® U(1)y- This allows 
for the possibility that our effective theory is UV completed by a theory which has already 
undergone electroweak symmetry breaking, so that operators such as the d = 5 Higgs portal 
[25]-[29] operator xxY^V^ (where VV can be Z°Z° or V^ + VI^ _ ) may be considered here. 
Throughout this work we present bounds assuming Dirac x> as corresponding bounds for 
Majorana x differ only by the appropriate symmetry factors or vanish identically, as in the 
case of operators containing xi^X-> or XP^X- 

Since we are not model-building UV complete dark sectors, but rather working in an 
effective operator approximation our assumptions about the electroweak charges and repre- 

5 In this document we will use the symbol \ to refer to both the DM particle itself and to the multiplet 
containing the DM particle. When it is not clear from the context we further clarify using x° and x ± - 
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sentations of the \ multiplet have little qualitative impact on our end results. Nevertheless, it 
is conceptually helpful to discuss this here in some detail. DM must be the electrically neutral 
and stable component of some multiplet x an d may be either a SM singlet or non-singlet, ex- 
cept in the case of operators such as the "magnetic moment," X^fx^ GF ? which obviously 
require a x transforming under SU(2) L ® U{1)y- For multiplets having both charged and 
neutral components the charged states are generically heavier than the neutral component 
by 0(100 MeV — 1 GeV) due to loop corrections [30j[31j, although a full description of the 
UV theory would be necessary to actually calculate this splitting and to assure that x° is the 
lightest new particle. Additional discrete symmetries (e.g., R-parity, T-parity, KK-parity, 
etc.) are typically posited to prevent x° from decaying into lighter SM states. Non-singlet 
X's may be in either chiral or non-chiral representations of SU(2) L ® £7(l)y. DM in chiral 
representations must get its mass, like the SM fermions, via the Higgs mechanism and thus 
has a renormalizable coupling, to the Higgs. Such a coupling is not considered in the 

current analysis and could significantly modify our results (although models with additional 
chiral matter are highly constrained by precision observables). Non-chiral may be further 
classified as being in real or complex representations of SU(2) L ® U{l)y (see [32]). For real 
representations T 3 = Y = 0, implying that the leading xxYV interactions are t-channel x ± 
exchanges built out of vertices X°X ± W T (the canonical example being the SUSY wino). For 
complex representations x has 7^0 and thus couples to Z° at tree-level, a coupling which 
badly violates DM direct detection bounds unless it is highly suppressed. This suppression 
is usually accomplished with two or more such multiplets, adding operators which mix gauge 
eigenstates with differing hypercharge such that the lightest neutral component hardly cou- 
ples to the Z° [32j (the canonical example being the SUSY higgsino). In the current work 
we try to remain agnostic as to the particular details of the dark sector and seek to present 
our bounds so that they can be interpreted in the context of any particular UV theory. 

While the EFT formalism offers simplicity and model-independence, its disadvantage, 
relative to the study of fully-defined models, is that it is not generally applicable. The 
crux of this issue is that when experiments probe our contact operators with energies 
such that the SM and DM particle are no longer the only light degrees of freedom (i.e., 
~ A for operators with dimensionful coefficient A~ n ) we expect our EFT description 
to break down and the details of the UV theory to become important. For direct detection 
experiments looking for DM scattering, in which characteristic energies are ~ O(10 keV), 
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the effective theory should essentially always hold (as evidenced by the typical x° ~ X ± mass 
splittings discussed above). For indirect detection experiments that look for the annihilation 
products of highly non-relativistic (v ~ 10 -3 ) halo DM, in which characteristic energies are 
~ 2m x , the effective theory should hold as long as the mass of particles that mediate 
these interactions is M J> 2m x , as we have basically assumed in stating that the SM and x 
are the only light degrees of freedom. One caveat to this is in cases where the masses of new 
particles which mediate the interaction are M < e 2m x ( "light mediators" ) . Although the EFT 
formalism may seem inappropriate in such a scenario, some studies have investigated taking 
this approach [10J [13J [15]. For collider experiments the efficacy of the operator formalism is 
not as automatic. In probing these operators via weak boson fusion at the LHC we expect 
energies on the order of the partonic CM energy ^ ~ C?(TeV). Strictly speaking 
then, bounds that are set in WBF searches must translate to bounds on mediator masses in 
excess of the average partonic CM energy in order to be interesting in this context. 

It is interesting to try and determine which portions of the A vs. m x plane cannot have 
a perturbative UV completion. This has been done in previous works, for example those 
studying XXQQ interactions, by identifying the dimensionful coefficient of these operators 
with the couplings and mass of an imagined heavy mediator exchange, 1/A 2 ~ gigz/M 2 , 
imposing M j> 2m x and g\g\ <; (47r) 2 , thus defining a perturbativity boundary in the A vs. 
m x plane. The situation is somewhat more complicated in the case of xxYV interactions. 
Consider, for example, the Higgs portal operator xx^^? which is a d = 5 operator with 
coefficient 1/A. If one imagines that the corresponding UV completion involves an s-channel 
exchange of some new heavy scalar particle S then the new vertex SVV in the UV theory 
is dimensionful, having a coefficient gsvs that may be associated to the electroweak vev 
v ~ 246 GeV, but could just as well involve dimensionful numbers generated via the details 
of the dark sector physics. Similarly, the coefficient of the d = 7 operator XXY^^^v is 
ambiguously connected to the parameters of the underlying UV physics, as such an operator 
typically describes a process that happens at loop level, having cross-sections which are 
complicated functions of the couplings and masses of particles in the loop (e.g., XiXi ~~ ^ 77 
in the MSSM). Given such complications we omit further discussion of UV perturbativity in 
this work. In any given UV theory this boundary can be straightforwardly computed and 
compared to the exclusion limits that will be presented below. 
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Unitarity offers a UV- insensitive criterion for determining the efficacy of our EFT 
description of WBF events at the LHC [33]. As is familiar from, e.g., pion scattering, 
amplitudes in a low-energy effective description may violate S-matrix unitarity as the energy 
of these interactions is increased. Conversely, at any given energy of interest we may typically 
consider increasing A so that, above some value, our effective description would not violate 
unitarity for interactions at this energy. For A below this value the apparent violation of 
unitarity signals the failure of our EFT description and suggests that some modification 
would be necessary (i.e., a proliferation of operators or the detailed dynamics of a particular 
UV completion) to correctly describe the physics. For events at the LHC we have the further 
complication that the energy flowing through our operator is distributed according to pdfs 
and kinematics. The most stringent use of unitarity on our EFT description would be to 
consider the worst case scenario, imagining the full machine center-of-mass energy flowing 
through our operator in order to set our "unitarity bound" ir^jthe A vs. m x plane. Typical 
WBF events, however, would have much less energy flowing through our operator because 
of pdf supression. We could then simulate events, calculating the energy flow through our 
operator in each, and derive contours in the A vs. m x plane such that, say, 99% of events do 
not violate unitarity. The interpretation of this is that most of the events that we include in 
our calculation of the experimental bound are well-described by our EFT language, only a 
small fraction are not (in an actual experiment these events might result in the production 
of particles that mediate this interaction and thus may not even look like WBF events) so 
the bound we derive should be fairly accurate. In what follows we will display contours for 
which 99%, 90% and 50% of event^] appear not to violate unitarity in 8TeV and 14TeV 
WBF events. 

Finally, while we will focus entirely on DM interactions with gauge bosons in this 
work, it is sensible to ask how the inclusion of the Higgs boson would affect our results. 
With SM singlet or vector-like DM, as is our focus, there is no possibility for renormalizable 
interactions coupling the \ to the Higgs before electroweak symmetry breaking. However, all 
combinations xX-> Xl^X-> e tc, are SM singlets and can be joined to operators like \H\ 2 , H*t a H 
and H*D^H giving non-renormalizable interactions between the \ an d the Higgs. After 

6 Keep in mind that this is not an experimental bound on the value of the UV scale in our low-energy 
effective description of DM, but rather, a bound on the values of A for which we can consider using our EFT. 

7 We do this calculation semi-analytically (details later) so our WBF analysis cuts have not been applied 
to these "events." The numbers 99%, 90% and 50% are only strictly correct then for operators with 100% 
signal acceptance. 
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electroweak symmetry breaking these operators lead to, e.g., xX^i XXYh, XXY an d XXYV 
interactions, the former two of which we are not considering in focusing on only operators of 
the form xxY or XXYV - There are two main consequences of omitting 3-point and 4-point 
operators with Higgses: in leaving out xxh interactions we will not have spin-independent 
scattering via the Higgs at tree level in our EFT, and in leaving out xxYh interactions we will 
not be able to quote gamma-ray line bounds for the 7/1 final state. Operators of the form xxY 
or XXYV arising from Higgses in the Lagrangian do not change our analysis qualitatively, 
but result in re-interpretations of these results. Taking xx\H\ 2 for example, gives a xxh 
coupling that (after integrating out the Higgs) would simply return what we are calling 
the Higgs portal operator: xxViiV^' Can this scenario be discerned from one in which some 
other (heavier) scalar mediates this interaction? At LHC energies, strictly speaking, the light 
Higgs cannot be integrated out, as it will actually be produced on-shell, but (as is usually 
done in light mediator analyses) one could use the narrow width approximation in order to 
derive collider bounds on xxVyV^ m either case. In considering direct detection bounds, the 
coupling to the light higgs XX^ could give significantly larger spin-independent scattering 
rates at a given A as compared to an interaction mediated by a heavy particle that does not 
couple to fermions directly at that same A (i.e., xxYyV 11 * s the only interaction in the EFT). 
Despite these being conceptually distinct scenarios, they likely cannot be distinguished in 
direct detection experiments and bounds on each are simply different interpretations of the 
same data. Combinations like H^D^H and (H* D ^H)(H* H) also return xxY an d XXYV 
interactions that, in the end, yield operators that we are already considering, e.g., 

ip^D^ip ^ ^(0 v)* d» -igW a ^ -ig'YB^ 1 ~ Z^. 



2.1 List of Operators 

We now discuss the particular operators that we will be working with. Table [T] lists these 
operators, according to their naive scaling dimensions, along with some of their properties. 
Column- by-column in Table [TJ we have listed operator names, Lagrangian expressions, our 
choice of canonical normalization, available vertex topologies in our EFT description of the 
operator, allowed sub-processes and the scaling of the leading terms in the non-relativistic 
expansion of the analytic formulae describing DM annihilation. 
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Name 


Expression 


Norm. 


Vertices 


Sub-Procs. 


Ann. 




dim = 5: 














D5a 


/V/V fJL 


A" 1 


4pt 


ZZ,WW 


V 2 




D5b 


XH5XV a ^ 


A" 1 


4pt 


ZZ,WW 


1 




D5c 




A" 1 


3/4pt 


A,Z,WW 


1 




D5d 


xv»XxV aiM/ 


A" 1 


3/4pt 


A,Z,WW 


1 (VV), v 2 


(//) 


dim = 6: 














D6a 


X^t a D uX V a ^ 


A" 2 


3/4pt 


A,Z,WW 


1 




D6b 


X1^5t a D u xV a ^ 


A" 2 


3/4pt 


A,Z,WW 


1 (VV), v 2 


iff) 


dim = 7: 














D7a 




A" 3 


4pt 


AA,AZ,ZZ,WW 


V 2 




D7b 


XHsXV^V^ 


A" 3 


4pt 


AA,AZ,ZZ,WW 


1 




D7c 


XXV^u 


A" 3 


4pt 


AA,AZ,ZZ,WW 


v 2 




D7d 




A" 3 


4pt 


AA,AZ,ZZ,WW 


1 





Table 1: Table describing operators used in this work. A full description is given in the text. 



The d = 5 operators D5a-b could result from exchanges mediated by new heavy 
scalar or pseudoscalar bosons. These operators necessarily require spontaneous breaking of 
SU (2) L ®U (l)y down to U (1) em and consequently only the and V1^ + V1^ _ subprocesses 

are allowed for these operators. The d = 5 operators D5c-d are similar to dark magnetic and 
electric dipole moment^] (respectively), which have been the subject of much recent study 
[34]-[4T]. As these operators have one field strength they give rise to both 4-point xxYV 
and 3-point xxY contact interactions. Since we are talking about interactions involving two 
neutral DM particles, the V^ v in the D5c-d operators can be either W 3fIl/ or B^ v ', giving 
specifically the subprocesses: XXW + W~ , XX%° an d XX^°- The d = 6 operators D6a-b [42] 
could arise via exchange of new neutral vector bosons (e.g., a Z f ) and, since = W 3fiiy or 
B^ v ', give rise to the same 4-point and 3-point interactions as in the D5c-d case. Finally, The 
d — 7 operators D7a-d typically arise from 1-loop diagrams (e.g., the XiXi ~~ ^ 77 process in 
SUSY), and may occur through any of the V1^ + V1^ _ , Z°Z°, Z°j and 77 subprocesses. 

8 Note that, due to the identity 75 cr^ = (z/2)e /izya/3 cr Q!/ 3, the operator D5d is equivalent to an operator of 
the form Xl^yLvXW^ '• 
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The operator expressions listed in Table [j] should all be understood as the sum of the 
expression listed and its complex conjugate expression. We employ canonical normalizations 
listed in the "Norm." column of the table. 

In the "Vertices" column we distinguish operators that have only 4-point interactions, 
from operators that have both 3-point and 4-point interactions. This is distinction is very 
important, as we will describe in the following sections, for both the LHC WBF bounds and 
for direct detection bounds. For WBF, operators with both 3- and 4-point interactions tend 
to generate events which look more like the background W/Z + jj processes than other- 
wise, and thus are somewhat harder to constrain at the LHC. As concerns direct detection, 
operators with only 4-point interactions must scatter via higher-order processes and cur- 
rent bounds have been estimated to be somewhat far from current experimental sensitivities 
[32j[43j[21j, while operators with both 3- and 4-point interactions may scatter off of nuclei 
via tree-level exchanges of the Z° or the A . Since we know that dark matter is dark and 
approximately collisionless, the 3-point coupling XX^ \ an d resulting long-range interaction, 
is highly constrained by an array of measurements [2] that are sensitive to both DM-SM 
interactions and to DM-DM self-interactions. The direct detection phenomenology of the 
XX A scenario diverges from that which is usually studied, as the recoil spectrum derived 
for such a long-range interaction is distinct from the recoil spectra derived for the typically 
assumed XX^N (SI) and Xl^lhX^l^l^N (SD) interactions, and has been the subject of 
recent developments [45j [46j [47] . Given that such a scenario is so tightly constrained we will 
suppose that the XX^° interaction is negligible in setting our bounds, e.g., by tuning La- 
grangian terms involving the W* v and B^ v fields to cancel this vertex. The only opportunity 

for appreciable scattering rates is thus the XX^° vertex. For our operators, however, these 
vertices always involve derivatives and are thus momentum suppressed and always result 
in negligible scattering. For m x <> mz/2 the XX%° vertex will allow the decay Z° — » xx> 
contributing to the invisible width of the Z°, which is measured to be r^ inv . <J 2MeV [48]. 
The constraint from this bound will be shown in the figures that follow. 

The entries in the "Sub-Procs." column refer to the allowed combinations of gauge 
bosons (W+Vl^ - , Z°Z°, 7Z and 77 - u subprocesses") that can arise from the generic and 
in each operator. Available subprocesses are determined simply by U(1)em invariance, 
as vertices need to conserve electric charge and photons are only allowed to arise from the 
field strength F^ v '. Any UV complete theory will specify the exact weighted combination of 
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subprocesses present in the EFT, but here we will need to make some assumption about these 
combinations in order to proceed. As is discussed in the following sections this assumption is 
more subtle in the context of WBF, where the gauge bosons are intermediate state particles 
and different subprocess combinations correspond to different coherent sums. One may also 
expect that custodial symmetry [49j [48j is an important constraint on the allowed subprocess 
combinations. We will not consider this in further detail, except to say that corrections (e.g., 
to gauge boson masses) are suppressed both by the appropriate power of the high scale A and 
also by an additional loop factor, so that combinations which deviate from the custodial limit 
are not obviously ruled out. In the end it is likely necessary that corrections to electroweak 
precision observables need to be calculated in the complete UV theory. 

The column labeled "Ann." lists the leading order terms in the non-relativistic ex- 
pansion of the analytic calculation of DM annihilation. For operators D5d and D6b this is 
quoted in terms of annihilations to either the VV or // final states, as they have different 
leading order terms in their expansions. We will obviously expect that operators with lead- 
ing order terms ~ v 2 will be significantly more difficult to constrain via indirect detection 
searches than those with leading order terms ~ 0(1). 

Finally, while not shown on the table, we note that the various operators have various 
P and CP properties. While parity violation in the dark sector is unconstrained (and possibly 
motivated by the observed parity violation in the SM), CP violation in the dark sector is 
possibly constrained via the CP violation that could be induced in the SM via higher-order 
interactions with DM in loops. Such constraints are expected to be highly model-dependent 
and are thus beyond the scope of this work. 

3 Bounds from Weak Boson Fusion at the LHC 

Weak boson fusion (WBF) processes have been widely studied as a means of enhancing LHC 
searches for the Higgs boson and such searches have been shown to be an effective strategy 
for discovering invisibly decaying Higgs boson^J Here we will use this same type of analysis, 
not for studying the Higgs, but for setting bounds on our contact operators. One could also 
imagine using searches for a single electroweak vector boson recoiling off of missing transverse 

9 with cross-sections, after cuts, of order of ~ 100 fb at the 14 TeV LHC 
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momentum, as has been discussed in the context of searches for an invisibly decaying Higgs 
in, e.g. ^ [50J. Here we focus only on WBF, as the enhanced production through longitudinally 
polarized gauge bosons has the potential for further reach than otherwise. For operators that 
favor transversely polarized gauge bosons, however, we expect that such a mono-V" searches 
would complement the results derived here. WBF signal events are characterized by the 
presence of two very energetic and well-separated forward/backward jets, as well as large 
missing transverse momentum. Here we derive bounds on our operators from WBF searches 
analyses of (current) 8TeV 25 fb -1 and (future) 14TeV 100 fb -1 LHC data sets. 

This section is divided up into three parts: first we discussed the application of the 
contact operator approach in the context of WBF searches, next we discuss the details of our 
WBF analysis in detail, and finally we present the resulting bounds on our contact operators. 

3.1 Contact Operators and Weak Boson Fusion 

As mentioned in the introduction, employing our EFT description in the collider environ- 
ment is much subtler than an EFT description of direct and indirect dark matter searches. 
The first issue that we need to address is the treatment of subprocesses in our WBF anal- 
ysis. Generic UV completions of our EFT will result in non-trivial combinations of contact 
interactions connecting the DM particle x to the gauge boson pairs ty + Vl^ _ , 7Z and 

77 (subprocesses). For indirect detection searches different combinations of final state gauge 
bosons will simply result in a different combinations of exclusion limits from the various DM 
searches. In a WBF search however, the gauge bosons are intermediate state particles, and 
in principle the different subprocesses add coherently. Here we have dealt with this difficulty 
as follows. 

In a particular UV theory one can explicitly calculate the coefficients that relate these 
different subprocesses, which we denote schematically as: 

otw+w~ ) a z°z°i a -yZ°i an d c^ 77 . (1) 

If only one of these a's is non-zero then we can obtain the exact bound by doing a col- 
lider simulation and calculating a cross-section for that particular subprocess operator in 
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isolation 
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°w+w-(s), vz°z°(s), cr 7 zo(^), ora 77 (s), (2) 

What we will do for this analysis is to simply present limits based on these individual 
subprocesses, and to estimate the limit on an operator involving a sum of sub-processes by 
the weighted incoherent sum of these limits: 

atat = ^o^(s) where > i £ {W+W~, Z°Z°, 7 Z°, 77}. (3) 



Of course, there should in principle be constructive or destructive interference amongst the 
amplitudes for different subprocesses, resulting in larger or smaller total cross-sections. It 
is clear that the procedure described above cannot account for this, however we find that 
the error incurred in employing this approximation is typically small compared to the other 
sources of error (e.g., systematic uncertainties) present in the estimation of the bound. 

To give a concrete example, suppose that the operator D7a is the relevant WBF 
operator resulting from integrating out the heavy fields in some particular UV theory. We 
then assume that the SU(2) L and U{l)y fields W a ^ v and B^ v appear with relative weights 
a and b. In this case the momentum space Feynman rules ("F.R." below) derived in the 
mass-eigenstate basis are: 

1 tp d An, 

^( ac w + j^i/i 2 (Z°Z°) 

X w , (77) 
(7Z ) 



A 3 




4(as 2 w + 


K) 


A 3 




Ac w s w (a 


-6) 


A 3 



where p 1? p 2 ? Mi an d 1^2 are the momenta and Lorentz indices of gauge bosons, s w and c w are 
the sine and cosine of the Weinberg angle, and X Ml/i2 = p^p^ 1 — p\ • p2 r r] filfJj2 . Here we define 
the cVs by taking the squared prefactors from the above equation (we actually generate the 



10 Here the (s) is a reminder that these cr's are dependent on the machine center of mass energy and also 
on the cuts employed in the analysis. 



12 



cross-sections Oi using A = 1 TeV, as is explained in Section 3.3): 



ITeV 
A 



6 




otw+w- — IGa 



(4) 



16(as w + bcJ I — ^— 1 



6 



. ,.. 2 /lTeV\ 



6 



'77 



CLy^o = 16 (Cj 



"ID 



Thus we compute for this operator as: 



&tot = 



(5) 



= 16 (^^— J {( ft ) 2 %+^-( 5 ) + ( ac l + &4) 2 °zozo(s) 
+ (as* + 6c^)V 77 (s) + (c^^(a - 6))V 7Z o(s)}. 




We will explore the bounds obtained for the example given above in further detail in Section 
3. 3| (after we describe our numerical analysis in detail), where we will pay particular attention 
to the influence of interference (or lack thereof) and systematic uncertainties on these results. 



3.2 Selection for WBF at LHC: SM Background generation 



There are several classes of important backgrounds to our signal, which are illustrated in 
Figure [I] We have: (i) Drell-Yan processes with two colored particles radiated off of the 
initial state serving as the two leading jets, i.e., Zjj and Wjj where Z° —> vv and W —> 
Iv (where in the latter case the lepton is not identified), these are referred to as QCD 
Zjj or QCD Wjj backgrounds, (ii) processes involving the t-channel exchange of a weak 
gauge boson and further radiation of another EW gauge boson, these will be referred to as 
EW Zjj or EW Wjj and (Hi) pure QCD backgrounds, where the mis-measurement of jets 
leads to significant missing transverse momentum. Of these, the pure QCD component can 
be substantially suppressed by cutting on p™ ss , so we simply omit these events from our 
background simulation (this was also found for the 3j simulations performed in the work 



This analysis is carried out at the parton level, following the analysis presented in 
[5T] as closely as possible. An important difference between the two analyses is that, in 



ISU). 
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(a) (b) (c) 

Figure 1: Example classes of feynman diagrams for Zjj background, similar diagrams are 
found for Wjj background as well. 

the work [5T] . the authors used various calculational tools to simulate the various classes 
of events, while our analysis makes exclusive use of the Madgraph v. 5 [52j package for all 
SM backgrounds, as well as for the signal events arising from the operators described above. 
While these differences are reflected in the cross-section limits derived in the two analyses 
we observe that these discrepancies do not have a qualitative impact on setting limits on the 
high mass scale, A. 

Here we describe the most relevant cuts that have been applied in our analysis and 
discuss their effects on the various background and signal components. First, we impose 
cuts to select events with hard forward and backward jets that are widely separated in 
pseudorapidity: 

p J T >40Gey, 1^1 < 5.0, 
\Vji ~ Vj2\ > 4.4, r]ji • r) j2 < 0. (6) 

The origin of these cuts in the Higgs WBF analyses traces back to the fact that Higgses 
produced in weak boson fusion are produced predominantly by longitudinally polarized W 
bosons [^J which are radiated preferentially with p T ~ giving high rapidity leading 

jets. In contrast, the jets produced in the background events are much more centrally 
distributed. Importantly, this latter fact is seen to arise because of interference between the 
two diagrams Figs. [TJd and[j]3. As we will discuss in more detail later on, this has important 
implications for the signal rates from some of our operators as well. Next, we impose a cut 
on the missing transverse momentum of events, £>™ S5 , in order to remove the contributions 
from QCD 3j and the soft single z/s resulting from W decays, as these contributions fall off 

11 This is more easily seen in the effective- W-approximation [53], where WBF is treated as a 2-to-2 process 
initialized by two W's, which are treated as approximately on-shell partons inside of the proton. 
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j] Distribution for 10K events for (Higgs, ZEW) Aij Distribution for 10K events for (Higgs, ZEW) 




-4 -2 2 4 -10 -5 5 10 

rj At? 



Figure 2: We display the jet rapidity distributions of a higgs-like signal (red) and electroweak 
Zjj background (blue) for samples of 10,000 events each (left panel.) Right panel shows the 
Arj = rjji — rjj2 distribution, we see that a cut of Arj > 4.4 roughly contains only the tails of 
the background distribution while including the full peak of the higgs-like signal distribution. 
Other background components behave similarly. 



very quickly with energy above ^ 100 GeV (c./., Fig. 1 of [51J): 

p™ ss > 100 GeV. (7) 

We further impose a cut on the invariant mass of the two tagging jets to suppress the 
contributions from QCD Zjj/Wjj background events, whose radiated gluon jets are typically 
softer than those of the corresponding quark jets in EW Zjj/Wjj (and signal) events (c./., 
the steeply falling da/dMjj of these events in Fig. 2 of [51]): 

M jlJ2 > 1200 GeV. (8) 

We further note that, since there is color exchange in the t-channel of the QCD Zjj/Wjj 
processes, these events tend to result in higher jet activity in the central part of the detector 
as compared to the EW Zjj/Wjj and signal events. To account for this we follow [5]] [54] in 
simulating the effect of requiring a pt > 20 GeV veto on jets in the central region by simply 
applying survival probabilities of 0.28 and 0.82 to QCD Zjj/Wjj processes and EW Zjj/Wjj 
processes, respectively, rather than actually applying the cut on an event-by-event basis. To 
remove a majority of Wjj background events, we veto events in which we can confidently 
identify the lepton according to the criteria: \rji\ < 2.5 and pri > 5, 10 and 20 GeV for e, /i 
and r leptons, respectively. Finally, we apply a cut on the azimuthal interval, A0, between 
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the two tagging jets, which is especially helpful in discriminating events according to their 
Lorentz tensor structure. We apply the cut 

A 0= \cj) 3l -cj) j2 \ < 1. (9) 

It has been observed [51j[55j[56j that such a cut favors Higgs-type contact interactions, 
whose jets tend to be relatively close in azimuthal angle, relative to the QCD Zjj/Wjj and 
EW Zjj/Wjj backgrounds whose jets tend to be more back-to-back in the azimuthal plane. 
In Figure [3] we give Acj) distributions for the various background components in this analysis, 
as well as for the various contact interactions used in this analysis. Since we are interested 
in signals coming from contact interactions with a variety of different Lorentz structures, 
ranging from the Higgs-like operators of the form D5a-b to Z'-like operators of the form 
D6a-b, we don't always benefit from the discriminating power of this cut. Nevertheless, 
we choose to apply this cut in all cases in order to facilitate comparison with the existing 
literature [51] . In any case, the effect of this cut on the resulting bounds is found to be 
relatively small. 

For comparison purposes, Table [2] shows the 14TeV cross-sections for the various 
background components generated in our study, along with those found in [5]] . We observe 
that our simulated QCD Zjj/Wjj rates are about 70-80% of the size of those found in [51] . 
while the EW Zjj/Wjj rates from the two studies match well. The authors of [5]] argue, using 
partially data-driven background estimates, that one can achieve a combined systematic 
uncertainty on the background in this analysis of 3.0(1.2) % for 10(100) fb~ x at a 14TeV 
LHC. These estimates seem to be optimistic however, and we opt to instead quote 95% C.L. 
upper limits on total signal cross section, a tot) as a function of the systematic uncertainty on 
the backgrounds in Figure [4j In this way the reader can estimate the power of such searches 
under any plausible assumption of what the systematic uncertainty should be. 

Results for 8TeV 25 fb~ l WBF analyses are derived similarly to the 14TeV case, 
using exactly the same cuts as in the 14TeV case. 
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Figure 3: We display azimuthal interval (A0) distributions for the background (left panel) 
and signals (right panel) used in this study. The background is split into components: QCD 
Zjj (black), QCD Wjj (purple), EW Zjj (green) and EW Wjj (orange). In the right panel 
curves are unit-normalized, the signal curves being described by red, green and blue curves, 
along with the background EW Zjj shape (black-dashed) for comparison. The operators 
D5a-D5d and D6a-D6b all correspond approximately to the red curve, while the operators 
D7a-D7b correspond to the green curve and the operators D7c-D7d correspond to the blue 

curve. 



a (fb) 


QCD Zjj 


QCD Wjj 


EW Zjj 


EW Wjj 


Total 


[51J 


Here 


[51J 


Here 


[51j 


Here 


[51j 


Here 


[51J 


Here 


Eq 


s. (6 8) 


1254 


1055 


1284 


906 


151 


148 


101 


85 


2790 


2194 


Eqs. (6 


9) + C.J.V. 


71.8 


56.6 


70.2 


47.3 


14.8 


14.6 


9.9 


8.2 


167 


127 



Table 2: Total cross section (in fb) for various background components after applying par- 
ticular sets of cuts. The results found in the work [5l] are presented along with those found 
in this work for comparison. Here "C.J.V." refers to application of survival probabilities for 
the central jet veto, as described in [5lj . 
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95% CL Upper Limit on cr tot 95% CL Upper Limit on cr tot 




5 10 15 20 5 10 15 20 

% BG Sys. Uncertainty % BG Sys. Uncertainty 



Figure 4: Here we display the 95% CL upper limit on the total signal cross-section, a tot 
(fb) as a function of an assumed level of systematic uncertainty in the background. In the 
left panel we compare the results found in this work and those that would be found using 
the backgrounds derived in the work [51] for 14TeV WBF analyses with lOOfb -1 (solid) or 
10fb -1 (dashed) data sets. In the right panel we compare our limits for a 14TeV lOOfb -1 
analysis with those for a 8TeV 25fb _1 analysis. 



3.3 Selection for WBF at LHC: Signal Subprocess Cross-Section 

Signal cross-sections are generated using the same procedure as that described in Section 



3.2 for generating SM backgrounds. We use FeynRules [57f^| in concert with Madgraph v5 
[52] to simulate signal events from our effective contact interactions. We have generated 
all signal subprocess cross-sections with A = 1 TeV. The signal cross-sections are oc 1/A^ n 
with n=l, 2 and 3 for contact operators of dimension d = 5, 6 and 7, respectively, so we 
can obtain results for A s ^ 1 TeV by the appropriate scaling. Subprocess cross-sections (for 
A = 1 TeV) and 95% CL lower limits on A are displayed in Figures [5]{9| 

In cases where an operator is associated with multiple subprocesses we have to choose 
a definite weighting for the subprocesses in order to plot limits on A. The combinations that 



12 FeynRules is a Mathematica package that allows for the calculation of momentum space Feynman rules 
for quite generic models of new physics (specified at the Lagrangian level). Here we have been able to 
interface FeynRules and Madgraph v5 using the Universal FeynRules Output (UFO) language. 
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were chosen to produce Figs. [5}(9]are: 



1 f Z^Z \ 
D5a: -%X [—^ + W+»W~ + h.c. j (10) 

1 / Z^Z \ 

D5b : jXH5X l—^ + W + »W~ + h.c.) 

D5c : ^ (x°^xW^ + ^X<W^) 

D5d : f (x<V* 3 X^ + ~^X^XB^ 

D6a : ^ (xT^X^ + *7,A,X^) 

D6b : ^ (x757m^^X^ 3 ^ + ^X757mA,X^) 

D7a : ^XXW^W^ 
D7b : i_^ 75X ^^ 

D7c : ^W^W^ 
D7d : ^WsxW^W* 



The combinations shown for operators D5c-d and D6a-b are fixed by the requirement that 
the xxA° vertex vanishes but all of the other expressions have been chosen arbitrarily for 
simplicity. In general one can weight the cross-section curves given in Figs. [5]|9] and sum 
to get a total cross-section according to the details of their own UV complete theory as 
described in Section 13.11 

The general features of the curves in Figs. [5j|9]are easy to understand: sensitivity is 
constant for all DM masses m x <C and falls dramatically at m x ~ 1 TeV ~ y/S. There is 
little difference in pairs of operators that differ only by a 75, except where the phase space of 
the DM particles becomes important near m x ~ The differences in overall reach between 
the various operators arise not only from the basic fact that they have different naive scaling 
dimensions, but also because of kinematic differences that have an impact on the effectiveness 
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of our WBF cuts. The operators D5c-d and D6a-b are all operators for which we have both 
3-point and 4-point contact interactions. The 3-point vertices are problematic here, as they 
result in diagrams which are kinematically similar to the SM backgrounds (Fig. [T]). In these 
cases we expect our WBF search to perform significantly worse as our cuts are not as able 
to distinguish signal and background. In Figs. [5}{9] we also show the regions constrainecp*] by 
the measurement of the invisible width of the Z°. 

Although our d = 7 operators do not suffer from the above issue, their WBF signal 
events are also somewhat difficult to separate from the backgrounds. This happens because 
the D7 operators are not predominantly produced from longitudinally polarized W bosons, 
which are both the dominant source of electroweak gauge bosons in the beam protons and 
are preferentially supplied with low transverse momentum. Thus the signal event rates are 
somewhat lower and the leading jets from WBF production of the D7a-d operators tend 
to be more central, suppressing signal efficiencies for the leading jet cuts. For illustration 



leading jet rapidity distributions are shown for many of our operators in Figure 10 



In Figures [5j|9] we also display curves describing the unitarity of our effective operator 
description. We include curves for which approximately 99%, 90%, or 50% of simulated 
events appear not to violate unitarity (red and pink curves are appropriate for 14 TeV and 
8TeV events, respectively). Roughly these curves separate regions above which our EFT is 
an excellent (99%), good (90%) or poor (50%) description as far as unitarity is concerned. 
These curves were calculated semi-analytically, combining the closed-form longitudinally 
polarized amplitudes for each operator (as calculated using FeynArts 3.4 [58] ) with the 
approximate analytical parton luminosity for longitudinally-polarized W-bosons derived in 
[53] and CTEQ5M parton distribution functions. Although the unitarity curves for operators 
related simply by a 75 are somewhat different we show only the curves associated to the "non- 
75" operators in Figures [M8] in order to reduce clutter (curves for the operators with 75 's 



are displayed later on in Figs. 22 25). We see that the Higgs Portal operators D5a,b are the 
most challenged in this regard, with the WBF excluded regions reaching above the "good" or 
"excellent" curves only for very light m x <J 50 GeV. The dipole moment and vector exchange 
operators D5c,d appear to be good descriptions up to m x of several hundred GeV and the 
d = 7 operators have no problems with unitarity up until the kinematic reach of the searches 



13 Since the requirement that our DM is not milli- charged fixes the relationship between the coefficients of 
SU (2)l and hypercharge gauge bosons in our EFT description one has no way of getting around this invisible 
width constraint unless somehow the very tight constraints [44 on milli-charged DM can be avoided. 
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m x ~ 1 TeV. 
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Figure 5: We display subprocess cross-sections (with A = 1 TeV) for the d = 5 "Higgs-portaF 
operators D5a (solid) and D5b (dashed). W+W~ and Z°Z° subprocesses are described by 
blue and red curves, respectively. Results for 14 TeV and 8 TeV scenarios are shown in the 
left and center panels, respectively. In the right panel we display curves representing 95% 
CL lower limits on A due to the subprocess combinations described for these operators in 
Eqn. 10, Bounds are set assuming either 14 TeV 100 fb _1 (black) or 8 TeV 25 fb _1 (grey) 
WBF analyses. In both cases limits are computed assuming 5% systematic uncertainty on 
the background. Red and pink curves describing the unitarity of the EFT description are 
also included, as explained in the text. 



Now let us follow up on our earlier discussion (Section 3.1) of errors and uncertainties 
in our WBF analysis, employing as an example the operator D7a. We would like to illustrate 
the effect that the various sources of error/uncertainty that we have been discussing, e.g., the 
incoherent sum of subprocesses, systematic uncertainties on the background and differences 
between the backgrounds as calculated here and elsewhere, in the context of this example. 
The panels of Figure [TT] address each of these effects in turn. 



In the upper-left panel of Fig. [II] we display subprocess cross-sections with weighting 
(eg j OLiGi, with as in Eq. ([5])), taking for simplicity a = 1 and 6 = 0. We also display the 
total cross sections resulting from coherent and incoherent summation of the subprocesses, 
o coherent and ^incoherent- I n the upper-right panel we convert these total cross-sections into 
95% C.L. lower limits on A, ^coherent and incoherent- We observe that there is constructive 
interference between the different subprocesses in this example, resulting in a fractional error 
of « 20% in using the approximation a tot « ^incoherent- We note, however, that the fractional 
error induced in using the bound incoherent instead of Coherent is related to the fractional 
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Figure 6: Similar to previous figure but now for the d = 5 "moment" operators D5c (solid) 
and D5d (dashed). For these operators the only relevant subprocess is DM production 
through W + W~ and can occur via either 4-point or 3-point interactions. In the right panel 
we plot the limit from the Z° invisible decay width (purple) along with the results for 14 TeV 
100 fb" 1 (black) or 8 TeV 25 fb" 1 (grey) WBF analyses. Red and pink curves describing the 
unitarity of the EFT description are also included, as explained in the text. 
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Figure 7: Same as in the previous figure but now for the d = 6 operators D6a (solid) and 
D6b (dashed). 



22 



D7a (Solid) and D7b (Dashed) D7a (Solid) and D7b (Dashed) D7a (Solid) and D7b (Dashed) 




m x (GeV) m x (GeV) m x (GeV) 

Figure 8: We display subprocess cross-sections (with A = 1 TeV) for the d = 7 operators 
D7a (solid) and D7b (dashed). The subprocesses V^ + VI^ _ , Z°Z°, 77 and jZ° may all result 
from this class of operators and are represented here by blue, red, green and purple curves, 
respectively. Results for 14 TeV and 8 TeV scenarios are shown in the left and center panels, 
respectively. In the right panel we display curves representing 95% CL lower limits on A 
due to the subprocess combinations described for these operators in Eqn. [TOj Bounds are 
set assuming: either 14 TeV 100 fb" 1 (black) or 8 TeV 25 fb" 1 (grey) WBF analyses. In both 
cases limits are computed assuming 5% systematic uncertainty on the background. Red and 
pink curves describing the unitarity of the EFT description are also included, as explained 
in the text. 
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Figure 9: Similar to the previous figure but now for operators D7c (solid) and D7d (dashed). 
Unitarity curves were not calculated for these operators but are expected to be similar to 
those found for the operators D7a,b. 
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Figure 10: We display the jet rapidity distributions of various operators for 10,000 events, 
i.e. they are not normalized to their respective cross-sections. 



error in the total cross-section as 



6 A _ 1 5a 



(again, n=l, 2 and 3 for contact operators of dimension d = 5, 6 and 7, respectively) so that, 
e.g., a 20% error due to neglecting interference in a tot will induce an error in A of only 5%. 



This is reflected in the upper-right panel of figure 11 



In the lower-left panel of Fig. [TT] we display 95% C.L. lower limits, in all cases using 
^incoherent) for four different values of assumed systematic uncertainty on the background 
(0%, 5%, 10% and 20%). The variation in these curves is significantly larger than that seen 
in the other panels. 



In the lower-right panel of Fig. 11 we display 95% C.L. lower limits (in all cases 
using ^incoherent and ignoring the systematic uncertainty) employing either the background 
rates found in this analysis or the background rates found in the work [51J. As mentioned 
previously, for a given assumed systematic uncertainty, the difference in bounds resulting 
from using the backgrounds derived in [51 J and in using the backgrounds derived here is 
relatively insignificant. Overall we see that systematic uncertainty is expected to be the 
dominant source of error in our analysis. 
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Figure 11: We display WBF bounds for 14TeV/100 fb _1 LHC analyses on the operator 
D7a. In the upper-left panel we display weighted subprocess cross-sections (as described in 
the figure) along with the coherent (black-solid) and incoherent (black-dashed) total signal 
cross-sections. In the upper-right panel we translate coherent and incoherent total cross- 
sections into bounds on the high scale A. In the lower-left panel we describe the effect 
of background systematic uncertainty on the A bound, showing 0%, 5%, 10% and 20% 
uncertainty, as denoted in the figure. In the lower-right panel we show the difference in 
bounds that can be placed using the backgrounds derived in this analysis (orange-solid) or 
in using the backgrounds derived in the work [5]] (orange-dashed). 
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4 Dark Matter Search Bounds 



Here we investigate bounds on our contact operators that can be derived from astrophysical 
data. As we discussed in Section [2j signals in experiments which probe our contact inter- 
action via non-relativistic DM scattering or annihilation are particularly well-modeled by 
an effective operator description, in possible contrast to WBF bounds derived in this con- 
text. Of course, DM signals from astrophysical DM distributions are also subject to many 
sources of uncertainty, e.g., in estimating average DM relic abundance [59J, local effects of 
DM substructure [6Q]-[64J, propagation of SM products of DM annihilation, etc., that are 
not present in collider searches. As we discussed earlier, we do not expect significant direct 
detection scattering rates from our operators as operators with only 4-point contact interac- 
tions must scatter through higher order processes and as operators that generate the XX% 
three-point vertex are momentum suppressed. Given this, our primary focus here will be 
to derive bounds that can be set from null searches for DM annihilation in our Milky Way 
(MW) DM halo. 

Earth-bound and satellite-born detectors search for the products of DM annihilations 
in the MW halo by measuring a variety of energetic particle spectra. Our operators produce 
a wide variety of SM final states (W + V1^ - , Z°Z°, 7Z and 77 for operators with 4-point 
interactions only and additional fermionic final states for 3-point interactions where DM 
annihilates through s-channel 7/^° ) and thus can be constrained by a variety of indirect 
detection experiments. Here we focus on three classes of search: 7-ray spectral limits from 
MW dwarf spheroidal galaxies, monochromatic 7-ray line searches in the MW halo, and 
measurements of antiproton cosmic-ray spectra. 

4.1 Indirect Detection Bounds 

Dwarf spheroidal galaxies are extremely DM dominated satellites of the MW which, having 
essentially zero intrinsic astrophysical 7-ray sources, are an excellent place to look for the 
continuum 7-ray spectra that accompanies DM annihilation into essentially all final states. 
For this bound we employ the 95% C.L. limits set by the Fermi-LAT collaboration using a 
combination of observations of ten MW dwarf-spheroidals [65j , as well as limits derived from 
the VERITAS collaboration's observations of the MW dwarf Segue-I [66]. The low-energy 
threshold of the satellite-born LAT instrument extends far below the ~ 100 GeV threshold 
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of the earth-bound VERITAS air Cerenkov telescope array so that the LAT sets the most 
stringent limits for DM masses below ~ 900 GeV. Above this, the large fiducial volume of 
VERITAS (the atmosphere above the array) results in the tightest constraints on heavier 
DM. As the results quoted in [65j[66j are described in terms of limits on cross-sections into 
specific SM final-state channels (WW, 66, rf and jijl) we do not have to model the signal 
7-ray spectra or the DM distributions in the dwarfs here. As has been observed [4], the 
continuum 7-ray spectra from annihilations to ZZ lead to essentially the same limits as 
that from annihilation to WW and that annihilations to all light quarks produce essentially 
the same limits as that from annihilation to 66. Given this, we sum the cross-sections for 
annihilation to the WW and ZZ final states and those for all light quarks in comparing 
to the experimental WW and 66 limit curves, respectively. In order to calculate the relic 
density, total annihilation cross-section and cross-sections for particular final state channels 
we use FeynRules 1.6.0 to calculate the Feynman rules for each operator and interface with 
MicrOMEGAs 2.2 [67J to calculate the DM observables. 



We use limits on 7-ray lines that were set by the Fermi-LAT collaboration [68], to 
bound operators that can annihilate directly into the 77 and jZ° and final states. We assume 
an NFW profile [69] for the MW DM halo. Results for the more conservative isothermal 
profile would be about 30-40% less constraining. Relating the DM mass to the 7-ray line 
energy, the LAT data provide constraints on 7 GeV < ttidm < 200 GeV DM annihilating to 
77 and on 60 GeV < tudm < 210 GeV DM annihilating to 7Z . 

We use the PAMELA collaboration's measurement [70] of the ratio of cosmic-ray 
antiprotons to protons to bound annihilations producing substantial hadronic matter. We 
again assume an NFW profile for the DM halo in calculating our antiproton signal rates and 



setting bounds. We use a modifiec^ version of the numerical package DarkSUSY 5.0.4 to 
calculate the signal antiproton injection spectra. This injection spectra is then propagated 
to obtain a local signal flux spectrum using a propagation model (galdef_50p_599278) that 
is supplied in the GALPROP v50.1p package [71] [72] [73] and is seen to be a good fit to a 



variety of astrophysical observations. We calculate the bounds from the PAMELA data 



15 



by calculating a x 2 , where we only include signal protons/antiprotons and dividing by the 
(well-measured) primary proton cosmic-ray spectrum, excluding regions at 95% confidence 



14 modified to use non-SUSY models. 

15 Taking the 17 highest energy bins, as the lower energy bins are affected by solar modulation. 
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(x 2 /16 > 1.724). As we are not adding any astrophysically produced "secondary antipro- 
tons" this exclusion is somewhat conservative, though calculations done with background 
added are seen to provide similar bounds. It should be noted that the numerical tables 
used in DarkSUSY to derive the injection spectrum were created by scanning DM masses 
in the range lOGeV — 5TeV, so one has less confidence in the numerical accuracy when 
extrapolating beyond this range (e.g., in looking at light dark matter scenarios) 16 , 

In using astrophysical experiments to set bounds on our operators we must make some 
assumption about the relic abundance of our DM. In all of the figures that follow we deter- 
mine excluded regions by assuming that the dark matter relic density is Qh^ = 0/i^ MAP ~ 
0.114 (i.e., that the x particle makes up all of DM). In the figures we include a (red-dashed) 
curve in the A vs. m x plane where the annihilation through the contact operators alone 
would yield a relic density, calculated assuming the usual thermal cosmological evolution, 
that matches the WMAP [76] value. In any fully specified model, however, the true relic 
density of the DM may be greater or less than this thermal value. For example, the relic 
density could be decreased by annihilating to dark sector states that later decay into the 
SM, or increased by non-thermal cosmological evolution [77j[78j[79j or co- annihilation with 
states that have relatively inefficient annihilation and are not included in our EFT picture 
[80j[81j[82j. The current relic density of \ cannot be larger than the WMAP value (or DM 
would overdose the universe), but may be smaller than the WMAP value, in which case lim- 
its coming from searches for annihilating DM would be substantially less-constraining than 
what we show in our figures (annihilation signals scale as (f]^/f]^ MAP ) 2 ). If x is imagined 
to be a Dirac fermion then the relative abundance of \ an d X is important for determining 
limits. In this case we are assuming that the current relic abundance of x is equal to the cur- 
rent relic abundance of x- If this is not the case, and either the x or X abundance dominates 
(as in asymmetric DM scenarios [83]-[88]). then any bounds from annihilation would vanish. 
To zeroeth-order, the reach of indirect detection limits on a given operator just depends on 
whether the operator is velocity suppressed in the non-relativistic limit. Such operators are 
scaled by a factor v 2 ~ 10 -6 compared to unsuppressed operators. A characterization of the 
leading order scaling of each operator can be found in the "Annihilation" column of Table 

16 Of course, one should also consider the range over which the underlying code that has generated these 
tables (PYTHIA V6 [74]) is accurate. The authors of DarkSUSY note, in particular, that there is something 
like a factor of 2 uncertainty in the antiproton yields at low energies due to the lack of reliable low-energy 
antiproton data with which to tune PYTHIA [75] . 
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□ 



4.1.1 Results 



Figures [12|f2T] describe relic density contours and indirect detection bounds for our list of 
operators in the A vs. m x plane. The red-dashed line in all figures denotes where the 
thermal relic density calculated for the given operator (in isolation) matches the cold dark 
matter density measured by WMAP. Black lines in the left panel of all figures denote lines 
of constant relic density and may be used to scale annihilation signal rates to those expected 
in any given UV complete DM model. The region which is excluded by one or more searches 
combined is shaded in grey in the left panel of all figures. Additional panels show the regions 
excluded by each particular search separately, e.g., dwarf bounds on annihilation to /x-pairs 
or antiproton bounds from the PAMELA pjp measurement (of course, not all searches are 
relevant for any given operator). The excluded regions shown in Figs. 



12 



21 



were calculated 



for the particular operator combinations quoted in Eq. [10} As the gauge bosons are now in 
the final state for annihilation, there is no issue of coherence effects and it is fairly easy to 
derive excluded regions for other subprocess combinations, given the results presented here. 
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Figure 12: Relic density (left panel) and annihilation cross-section (right panel) contours in 
the A vs. m x plane for the Higgs-portal operator D5a. Contour labels are values of log 10 (Q/z^) 



and of log 10 ((crL>) cm 3 s -1 ), respectively. This operator is heavily velocity suppressed, so that 
even under the assumption that Qh^ ~ 

much too low for any exclusion from indirect detection searches. 



^wmapi the current annihilation cross-section is 
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Indirect Combined Dwarf Bounds p Bounds 

4 -°f/6 1 4 -°f 1 4 -°f 




Log m [m x (GeV)] Log m [m x (GeV)] Log m [m x (GeV)] 



Figure 13: Indirect detection limits for the operator D5b in the A vs. m x plane. The left 
panel includes contours of constant relic density (as calculated assuming thermal cosmology), 
labeled with values of log 10 (fi/i^), with the contour having Qh^ « ^h^ VMAP represented by 
the red-dashed line. The region which is excluded by one or more indirect searches is shaded 
in grey. In the middle and right panels we describe regions that are excluded by individual 
indirect limits. 
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Figure 15: Same as in Figure 13 but for the operator D5d. 
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Figure 16: Same as in Figure 13 but for the operator D6a. 
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Figure 17: Same as in Figure 13 but for the operator D6b. 
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Figure 18: Same as in Figure 13 but for the operator D7a. 
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Figure 19: Same as in Figure 13 but for the operator D7b. 
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Figure 20: Same as in Figure 13 but for the operator D7c. 
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Figure 21: Same as in Figure 13 but for the operator D7d. 



Many of the features of Figs. [12pT1 are easy to understand. The basic reach of 
the excluded regions varies quite a lot among the different operators according to velocity 
suppression intuition. Operators D5a, D6b, D7a, D7d are all velocity suppressed to some 
extent (none of the plotted region is excluded for the Higgs portal operator D5a), and so are 
relatively difficult to exclude. Operators with 3-point couplings to the Z° feature a prominent 
peak in sensitivity near m x ~ (while the relic density contours simultaneously peak). 

The various searches have effective sensitivities within obvious bounds: annihilation to WW 
only occurs for m x > raw, line searches are only effective for DM masses that generate lines 



of energy between the Fermi-LAT analysis thresholds^ and antiproton bounds reach as far 
down in m x as there are channels that produce hadronic matter (WW, 66, rf). Looking 

at the dwarf bounds, we see that excluded regions for WW and 66 final states feature a 
discontinuity around m x ~ 1 TeV, where the bound transitions from the Fermi-LAT limit 
to the VERITAS limitB 

Recall that the excluded regions for all searches are calculated under the assumption 
that Qh x = Qti^MAP an d that annihilations in the current epoch occur only through our 
contact operators. We see that this scenario is excluded for operators D5c and D6a for 
m x & via the Fermi-LAT dwarf 66 limit, and for operators D7b-c for m x mw, 

via the Fermi-LAT 7 line limits. In the former case this is a model-independent statement 



17 ' E 1 = m x for the 77 final state and E 1 = (m x — m 2 z /Am x ) for 7Z. 

18 We also observe that there is no such discontinuity in the fifl and rf excluded regions. The reason for this 
is that the Fermi-LAT and VERITAS limits on stiff spectra (from /ip and rf annihilations) approximately 
match up at their overlap, whereas the limits on softer spectra (WW and bb) do not. 
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(upon requiring that \ is n °t milli- charged), and this exclusion complements the constraint 
rz,inv. ^ 2MeV in this region (as we will see in the next section). In the latter case this 
is a model-dependent statement, as models that give D7b-c operators with negligible 77 



annihilation are not excluded. In the case where excluded regions in Figs. [T2jj21]do not reach 
up to the red line, where a standard thermal cosmological calculation using our operator 
would give Qh^ = fi^ M ^P5 then the excluded scenarios have a larger current annihilation 
cross-section than what would give Qh^ = ^/^map m a thermal cosmology. Such scenarios 
would require either (i) non-thermal cosmological evolution or (ii) inclusion of (relatively 
inefficient) annihilation channels other than our operators to boost the relic density of our 
X up to the WMAP value. 



5 Combined Results 



In Figures 22 26 we combine regions that can be excluded at 95% C.L. by 8TeV 25/6 -1 
and 14TeV 100/6 -1 weak boson fusion searches at the LHC, regions that are excluded at 
95% C.L. by current indirect detection searches and regions that are excluded by the Z° 
invisible width constraint, in the A vs. m x plane. For most operators WBF searches reach 
A ^ 500 GeV — 2TeV for DM masses m x <> ITeV, though operators that give three-point 
vertices are much harder to constrain in WBF, constraining A > 100 GeV only for DM 
masses below a couple of hundred GeV. Indirect detection reach varies wildly depending 
on the operator (predominantly according to velocity suppression, or lack thereof) and it is 
worth noting that these searches provide the only constraints for m x J> 1 TeV. For applicable 
operators, regions excluded by the invisible width of the Z° provide the tightest constraints 
for m x < mz/2. 

It is interesting to compare excluded regions to the overlaid relic density contours 



shown in Figs. 22 26, and to discuss what this means for our assumptions that (i) the DM 
abundance matches the WMAP value flh x « 0.11, (ii) that current annihilations proceed 
dominantly through our contact operators. Recall that the red-dashed lines in our Figures 
is the curve A(m x ) such that the thermal relic density, calculated using our operator in 
isolation, gives ~ 0.11. For operators D5c, D5d, D6a and D6b, the Z° invisible width 
constraint excludes such a scenario for m x <J mz/2, implying that DM of this kind with 
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™>x ^ m z/2 necessarily requires additional operators or dark sector annihilations in order 
to avoid overdosing the universe. WBF searches could similarly exclude this scenario for 
operators D7a-d with m x <J 100 — 200 GeV. For heavier m x we can only exclude scenarios 
with relatively larger annihilation cross-sections (i.e., scenarios that would require non- 
thermal cosmology, or something similar, in order to supply Qh^ « 0.11). 



D5a: WBF 95% Exclusion D5b: WBF and Indirect 95% Exclusions 




Log 10 K ( GeV )] LogioK ( GeV )] 



Figure 22: Regions excluded by 14TeV LHC WBF searches (8TeV- blue region, 14TeV- 
light blue region) and current indirect detection searches (green region) are compared in the 
A vs. m x plane for operators D5a (left panel) and D5b (right panel). Relic density is shown 
using the same contours that were described in Figs. [T2pT| Unitarity curves as described 
before for Figs. [5j|8]are reproduced here as well. 

6 Discussion 

In this work we have investigated bounds on fermionic dark matter which interacts with 
the Standard Model primarily via electroweak gauge bosons, using a model-independent 
effective field theory description. In this picture we can easily compare the experimental reach 
of different classes of DM search experiments, here focusing on the only relevant searches 
(given our assumptions): Weak Boson Fusion events in 8TeV or 14TeV LHC data, 7-ray 
observations of Milky Way dwarf spheroidals, 7-ray line searches, cosmic-ray antiproton data 
and constraints from the measured invisible width of the Z°. 

The combined reach of these experiments is seen to probe the UV cutoff scale, A, 
up to weak-scale values (hundreds of GeV up to several TeV). Scenarios in which dark 
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D5c: WBF, Indirect and r Z i nv 95% Exclusions D5d: WBF, Indirect and r Z i nv 95% Exclusions 




Log 10 K (° eV )] LogioK (° eV )] 



Figure 23: Regions excluded by 14TeV LHC WBF searches (8TeV- blue region, 14TeV- 
light blue region) , current indirect detection searches (green region) and measurements of the 
invisible width of the Z° (purple region) are compared in the A vs. m x plane for operators 
D5c (left panel) and D5d (right panel). Relic density is shown using the same contours 
that were described in Figs. [T2pl| Unitarity curves as described before for Figs. [5j|8] are 
reproduced here as well. 
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D7c: WBF and Indirect 95% Exclusions 



D7d: WBF and Indirect 95% Exclusions 





Log 10 K ( GeV )] 



Log 10 K ( GeV )] 



Figure 26: Same as in Figure 22 but for operators D7c (left panel) and D7d (right panel). 
Unitarity curves were not calculated for these operators but are expected to be similar to 
those in Fig. [25 , 
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matter interacts with Standard Model particles primarily via electroweak gauge bosons are, 
of course, harder to constrain than scenarios in which the dominant interactions are with 
strongly interacting particles. Nevertheless, the finding that searches are probing values of 
A near the weak scale is interesting as these are natural values for the strength of such 
interactions in the WIMP dark matter paradigm. 

We have discussed the implications of these bounds for the possible cosmological 
evolution of such dark matter, finding that relatively light dark matter scenarios (m x <> 
mz/2 or m x <J 100 — 200 GeV, depending on the operator) necessarily require additional 
structure (additional important operators or a non-trivial dark sector) to avoid overdosing 
the universe. 
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